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Determination of Radiative Widths of Scalar Mesons 
from Experimental Results on 77 — > tttt 
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The scalar mesons in the 1 GeV region constitute the Higgs sector 
of the strong interactions. They are responsible for the masses of 
all light flavour hadrons. However, the composition of these scalar 
states is far from clear, despite decades of experimental effort. The 
two photon couplings of the /o's are a guide to their structure. Two 
photon results from Mark II, Crystal Ball and CELLO prompt a 
new Amplitude Analysis of 77 — > 7r + 7r~, 7r°7r° cross-sections. De- 
spite their currently limited angular coverage and lack of polarized 
photons, we use a methodology that provides the nearest one can 
presently achieve to a model-independent partial wave separation. 
We find two distinct classes of solutions. Both have very similar 
two photon couplings for the / (980) and /o(400 — 1200). Hope- 
fully these definitive results will be a spur to dynamical calculations 
that will bring us a better understanding of these important states. 
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1 Introduction 



Two photon processes are a remarkably useful tool for studying the structure of matter 
and determining the composition of hadrons 0. Photons clearly couple to charged objects 
and the observed cross-sections are directly related to these charges. Thus, for example, 
in the reaction 77 — ► mi, the shape of the integrated cross-sections perfectly illustrates 
this dynamics. At low energies, the photon sees the pion as a whole entity and couples to 
its charge. Consequently, the cross-section for 77 — > 7r + 7r _ is large at threshold, whereas 
the 77 — > 7r°7T cross-section is very small 0. When the energy increases, the shortening 
of its wavelength enables the photon to see the individual constituents of the pion, couples 
to their charges and causes them to resonate (see, for instance, 0). Both the charged and 
neutral cross-sections are then dominated by the Breit-Wigner peak corresponding to the 
/2(1270) resonance, with several underlying /o states. The coupling of each of these to 
77 is a measure of the charges of their constituents (to the fourth power) and so helps to 
build up a picture of the inner nature of these mesons. But how do we determine their 
77 couplings from experimental data ? 

In an ideal world, with complete information on all the possible angular correlations 
between the initial and final state directions and spins, we could decompose the cross- 
sections into components with definite sets of quantum numbers. From these, we could 
then unambiguously deduce the couplings to two photons of all the resonances with those 
quantum numbers, not only the ^(1270) but also the more complicated scalar resonances 
/o(980) and / (400 - 1200) (and at higher energies the / (1500) and /j(1710)) [|]. Un- 
fortunately, in the real world, experiments have only a limited angular coverage and the 
polarization of the initial state is not measured. This lack of information plays a crucial 
role in any analysis and affects the determination of the resonance couplings || . Thus one 
has to make assumptions of varying degree of rigour : for instance, in the ^(1270) region, 
assuming the cross-section is wholly 1 = .D-wave with helicity two ||. Estimates of 
the underlying 1 = scalar couplings are made from the small 7r°7r° cross-section at low 
energies, or the much larger 7r + 7r~ cross-section, etc J7|, These are mere guesses and 
the consequent results of doubtful certitude. 

The aim of the present treatment is to perform an Amplitude Analysis in as model- 
independent way as possible. To achieve this, we make up for our lack of experimental 
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information, firstly by analysing the charged (tt + tt~) and neutral (7r°7r°) channels at the 
same time, and secondly using severe theoretical constraints from Low's low energy the- 
orem, crossing, analyticity and unitarity JS], ||. The low energy theorem means that the 



amplitude for Compton scattering 771 — > 777 is specified precisely at threshold [ 10 1 . An- 
alyticity, together with the fact that the pion is so much lighter than any other hadron, 
means the Born amplitude, modified in a calculable way, dominates the 77 — > nir process 



in the near threshold regime |TT], ||. This provides the anchor on to which to hook our 
Amplitude Analysis, determining all the partial waves with both I = 0, 2 below 5 or 600 
MeV. Above this energy, unitarity adds further constraints. Each 77 — > tttt partial wave 
amplitude is related to the corresponding hadronic processes h — > titt. Below 1 GeV or so, 
h can only be tttc and the constraints are highly restrictive. Above 1 GeV, KK channels 
not only open, but open strongly. This means we must incorporate coupled channel uni- 
tarity and include inputs from tttt — > KK too. Above 1.1 GeV, the rjr] channel opens and 
above 1.4 GeV a series of multi-pion channels become increasingly important. Because the 
r\r\ threshold is relatively weak, and the KK channel the major source of inelasticity, we 
can reliably perform an Amplitude Analysis, incorporating just ittt and KK information 
up to 1.4 GeV or so. Above that energy, we would have to access information on 77 — ■> mr 
and 7T7T — > nix (with n > 4) too and the analysis becomes impracticable, at present. 

Since 1990, when the last amplitude analysis of 77 — > tttt was performed ||, new results on 
77 — > 7T + 7r~ from the CELLO collaboration [12], more detailed information on the scalar 
7T7T final state interactions and increased statistics in the Crystal Ball experiment |H| on 
77 — > 7r°7r° have become available. These provide the impetus for a new analysis. In 
addition, there has been much speculation about the nature of the scalar states in this 
region, their relation to the lightest qq multiplet, to multiquark states and to glueball 
candidates @§ (T^. In each case, their two photon width is a key parameter in 
this debate. Consequently, we need to put what we presently know about such widths on 
as a firm a foundation as possible. Hopefully, this will be a spur to two photon studies at 
CLEO, at LEP and at future 5-factories. Further, improvements in data should allow the 
77 widths of the /o(1500) and /j(1710) to be fixed too. With this as the long term aim, 
the present analysis will be able to limit the number of possible solutions previously found 
and obtain more stringent information particularly on the scalar sector below 1.4 GeV. 
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2 Formalism and parametrization 



The unpolarized cross-section for dipion production by two real photons is given by the 
contributions of two helicity amplitudes M ++ and M + _ (the subscripts label the helicities 
of the incoming photons) [p], [5] : 
da 1 



1 — Ami/ s 



dtt 128tt 2 s 

These two helicity amplitudes can be decomposed into partial waves as 



T 



M ++ (s,6,<P) = eVie^E^s) y 



JO I 



J>0 



M + .(s,9,4>) = eVT6^£F J2 (s) YM4>) ■ (3) 

J>2 

The partial waves Fj\ (A = 0, 2) are the quantities we want to determine. 

As explained in the Introduction, such an Amplitude Analysis is not possible without 
some theoretical input. The first key constraint is unitarity. This relates the process of 
two photons producing some specific hadronic final state to the hadronic production of 
these same final particles. Thus, for each amplitude with definite spin J, helicity A and 
isospin /, unitarity requires (as illustrated in Fig. 1 for h = tttt, KK) 

Im^j A ( 77 - vrvr) = ^Jxill - h)* Tj(h - tttt) , (4) 

h 

where the sum is over all hadronic intermediate states h that are kinematically allowed; 
Ph being the density of states for each such channel. We have dropped any dependence 
the hadronic partial wave amplitude Tj may have on helicity, as we shall, in practice, only 
be concerned with spinless final and intermediate states. Eq. (4) is, of course, linear in 
the two photon amplitudes, T . However, each hadronic amplitude Tj{h — ► tttt) satisfies 
the non-linear unitarity relation: 

Im Tj{h -»• tttt) = Ph Tj(h -> h')* Tj(ti -> tttt) . (5) 



This equation means that Eq. (4) is satisfied by [18 



fjxdi -»• kti) = a h JX r j( h 7r7r ) > ( 6 ) 



h 



where the a^' JA are functions of s, which are real above im threshold. Thus, unitarity 
relates the 77 — > tttt partial wave amplitudes to a sum over hadronic amplitudes with 



the same tttt final state weighted by the coupling, a, of 77 to each contributing hadronic 
channel. Clearly, this constraint is only useful when we have information on all of the 
accessible hadronic channels. This restricts the present analysis to two photon energies 
below 1.4 GeV, where tttt and KK channels are essentially all that are relevant, Fig. 1. 
Then 

.Fj A (77 -> tttt) = a 1 /" -> tttt) + a 1 /" Tj(KK -> tttt) . (7) 




K 



Y 71 Y 7C Y 7C 

Figure 1: The relation between the two photon amplitudes and those for hadronic 
reactions given by Eq. (7). 

T-invariance of the strong interactions means 

Tj(KK -> tttt) = 7/ (tttt -> KK) . 

The analytic properties of the J-j\{s) suggest the functions a(s) are smooth for s > Am\, 
aside from possible poles that can occur in well-defined situations that we will discuss in 
more detail below. Notice that the a n ,K(s) give the weight with which Tj^nii — > rnv) and 
Tj{nn — > KK), respectively, contribute to ^ 7 ja(77 ^n). The a^s) and ax{s) will be 
determined by fitting the experimental data on 77 — > 7T7t, as we describe in Sect. 4. 

In Eq. (7) the hadronic amplitudes are independent of the photon helicity A, since the 
channels involve only spinless pions and kaons. Below 1 GeV, where the KK channel 
switches off, Eq. (6) can be expressed even more simply as 

•^a(77 - tttt) = a ! / X Tj(irn - tttt) (8) 

where a is a real function for 4m^ < s < Am 2 K . Eq. (7) and Eq. (8) are, of course, 
consistent, since for Ami < s < Am 2 K : 

r/(7T7T -> M) OC T/(7T7T -> 7T7t) (9) 

with a real function of proportionality. The use of the representation Eq. (7), throughout 
the region we consider, will allow us to track through the important KK threshold region. 
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We stop at 1.4 GeV, since there multi-pion channels (as well as 7777) become increasingly 
important and the unitarity constraint, Eq. (6), more complicated and impossible to 
implement without detailed partial wave information on tttt — ► An, Qn, etc. 

Though unitarity imposes Eq. (7), for each spin and isospin, in practice the 1 = 2 am- 
plitudes are simpler, as a result of the final state interactions being weaker and the KK 
channel not being accessible. Consequently, the representation, Eq. (7), will only be used 
for the I = S and .D-waves. Let us deal with these in turn: 

• 1 = S— wave: The 77 partial wave amplitude JFq will be parametrized in terms of 
the two real coupling functions a®' 00 and a^ 00 (denoted by a®, a° K as a shorthand) 
and the hadronic S'-wave amplitudes Ts(tttc — > ttti) and Ts(ttti — > KK). The 
input for the hadronic amplitudes is based on a modification (and extension) of 
the .fT-matrix parametrization of AMP |T8|]. Briefly, the T-matrix is related to the 
K-matrix by 

K 

= T^K < 10 > 

where p is the diagonal phase-space matrix. In the case in which two channels are 
considered, we have 

where 



» = I/I - , (12) 



1 / 4mi± 1 / 4mi , s 

P2 = ~\l ^ + ~\ 1 ^- (13) 



2 V s 2 

In this notation, the convention is 



1 <-> 7T7T , 2 <-> fCf^. 

Coupled channel unitarity is then fulfilled by the if-matrix being real for s > 4m^. 
It is then the i^-matrix elements, Kij that embody the hadronic information. For 
our 1 = 5*-wave, we take the K\ solution from ref. above approximately 1.1 
GeV, which is characterized by having only one pole in the K-matrix in the 1 GeV 
energy region. For energies up to about 1 GeV, we supply the strong interaction 
amplitudes Ts(tt7t — > ttti) and Tg(jnr — > KK) as given by the "Solution 2" obtained 
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Figure 2: Strong interaction inputs for the 1 = S-wave: here we show the 
moduli of the amplitudes T, which are obtained from the amplitudes T by 
removing the Adler zero factor, T(s) = T(s)/(s — s ). 



in a further refinement of the |18| fit, reported in ref. [BJ. It includes a larger 
set of experimental data, particularly in the KK threshold region, which favour a 
parametrization of the /o(980) that allows for two poles in the T-matrix on differ- 
ent sheets corresponding to this state. These solutions are smoothly joined from 
one energy region to the other. The moduli of the amplitudes Ts(ti7t — ► ttti) and 
Ts(7itt — > KK) determined in this way are shown in Fig. 1. In fact, we have two 
X-matrix parametrizations of these, called ReVAMPl and 2 []19"| , in which the 
X-matrix elements, in addition to having a single pole, have either 2nd or 3rd order 
polynomials, respectively. 

The relation of these hadronic amplitudes to that for 77 — ► ttti is through the 
coupling functions a n> K- Since these functions have only left hand cuts, they must 
be real for s > 4m^., but should be smooth along the right hand cut, once obvious 
dynamical structures are taken into account. Such obvious structures are that the 
coupling functions have poles wherever any element, or other sub-determinant, of 
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the T-matrix has a real zero. For instance, Ts(ttti — > 7T7t) and Ts(tt7t — > i^K) have 
Adler zeros below threshold (at s = So), which are not present in the 77 amplitude. 
(That is why they have been divided out in the amplitudes plotted in Fig. 2) [] 
Similarly, detT(s) can (and in fact does) vanish at a real value of s = si, where 
si < Am 2 K because of Eq. (9). Unless the a's have poles at this point, with related 
residues, this constraint on the hadronic amplitudes unnecessarily transmits to the 
77 amplitudes. To avoid these unnatural constraints, the functions are parametrized 
as follows: 

S — So S — Si 



a K \ s ) v 



<(s) = — T s (7r7r^7r7r;s = Sl ). (14) 

S — So s — Si 

The fit then determines and a^- as polynomials in s and the constant param- 
eter v. It is important to stress that the poles in Eq. (14) do not appear in the 
77 amplitudes, of course. This representation, Eq. (7), for the 1 = S-wave has 
the flexibility needed to determine the details of the mechanism by which the scalar 
resonances, / 's, couple to two photons. 

1 = D— waves with A = 0, 2: here a simplification arises from the fact that the 
nil — > 7i7i and 7i7i — > KK amplitudes are proportional to each other, both being 
dominated by the /2(1270) resonance. Then the hadronic amplitude Td(tt7t — > tttt) 
is given by 

_ . . BR 771/ r(s) . . 

T D (nn _ ^ = _ m) _ s _ imfT{s) (15) 

where BR is the branching ratio of / 2 — > titi. Importantly, the width is energy 
dependent and given by 

r(s > = mA™ (16) 



with j3(s) = yl — Am^/s. The factor D 2 (s) incorporates threshold and barrier 
effects. Here, we take these to be given by duality shaping, with the scale set by the 



1 in general the tttt — * tttt and tttt —> KK amplitudes do not have their Adler zeros at exactly the 
same point s = sq. However, present data are not sensitive to small differences between their positions 
and so in AMP [pi, these were taken to be at the same point, for simplicity. 
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slope of the non-strange Regge trajectories (or equivalently the p-mass, m p ). Then 
, / 6m 2 6m* \ , ( 6m 2 6m* \ 



m 2 — 4m 2 ) (m 2 f — Am 2 ) 2 J \ (s — 4m^) (s — 4m^.) 2 y 

(17) 

The tttt — > KK amplitude is proportional to this and so (just as in Eq. (8)) we can 
write 

ftx = a°' 2 \s) T d (tttt -> tttt) , (18) 

even above KK threshold, where the a° ,2X (s) are again smooth real functions of 
energy to be determined by the fit to experimental data. To ensure the appropriate 
threshold behaviour for the 77 — > tttt amplitude, these D-wave coupling functions 
are parametrized by modifying the threshold factor, so 

oP> 2 \s) = -£= . (19) 



m-f and r tot are taken from the PDG Tables 0] and D 2 (s) as in Eq. (17). 

So far the formalism we have described would apply to any reaction, by which a non- 
strongly interacting initial state leads to tttt. We now turn to the particular features of 
the two photon reaction. 



For 77 — > 7T7T, Low's low energy theorem [10] imposes an important constraint, in which 



the hadron charge fixes the size of the cross-section. This is embodied in the one pion 
exchange Born amplitude. Though the theorem applies at the threshold for the Compton 
process 77? — > 77?, the Born term controls the 77 — > tttt amplitude in the whole low energy 



region [11], as discussed extensively in 0. It is this dominance of the Born term that 
means, unusually for a strong interaction final state, that the 1 = 2 channel is just as 
important as that with 1 = 0. It is the almost exact cancellation between those amplitudes 
with 7 = 0,2 that makes the 77 — > 7r 7r° cross-section small at low energies. 

While the Born amplitude controls the low energy process, it is of course modified by final 
state interactions. As already mentioned these affect the 7 = and 7 = 2 tttt channels 
quite differently. For the 7 = 2 final state, these interactions are weak and the Born 
amplitude is little changed, remaining predominantly real in all partial waves. In contrast, 
the 7 = final state interactions in S and D waves are strong (leading to resonance 
formation for instance), consequently even close to tttt threshold the Born amplitude is 
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modified. It is unitarity that allows these modifications to be reliably calculated up to 



600 or 700 MeV |TTj, g]. Consequently, the Born amplitude, with such modifications from 
final state interactions, provides a precise description of the partial wave amplitudes on 
to which we must connect our amplitude analysis. 

• All waves with spin J > 4: For these, final state interactions are negligible and 
so the 77 — > it-it amplitudes are set equal to the Born amplitude in the whole energy 
region up to 1.4 GeV. Thus 

Fjx (77 -> 7T+7T-) = B JX (77 -> 7T+7T-) , (20) 

and 

Tjx (77 - ttV) = B JX (77 - ttV) = . (21) 



I = 2 S and D waves: Here the 77 — > tttt amplitudes have modifications from final 
state interactions that can be calculated up to 1.4 GeV. This we do by expressing 
the amplitude essentially as a modulus times a phase factor as 



1=2, 

JX \ 



where flj X 2 (s) is the appropriate Omnes function p0| : 



exp 



ds ■ 



(22) 



(23) 



with 



ij^ 2 the phase of the corresponding 77 amplitudes Fjx 2 - Applying elastic 



unitarity relates this phase to the 1 = 2 spin J tttt phase shift. Then, using the data 
of Ref. [pi]j , the Q's are readily computed. The function P(s) in Eq. ( p2|) is then 
real for s > 4m 2 and can be calculated as follows. Dropping the /, A indices to keep 
the notation simple, consider the analytic function fj(s) defined by 

Bj(a)(Slj\ a ) - 1), 



fj(s) 



(24) 



where Bj(s) is the spin J Born amplitude. We now write a once subtracted disper- 
sion relation for fj{s)/(s — 4m 2 ) J / 2 



fj(s) = fj(0) + 



+ 



s(s - Ami) 3 1 2 



7T 



4m2 s'(s'-s)(s'-4m^/ 2 



s(s-4rr4) J / 2 y° hnJgj^gT^gQ - 1)] 

7T /-oo s'(s'-s)(s'-4m2)J/2 ' 



(25) 
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where we have taken into account that Bj(s) has only a cut on the left-hand side, 
so that 

Im fj(s) = Bj(s)lm^lj 1 (s) for s > Ami > 

(26) 

Im fj{s) = Im [Bj(s) (fij^s) - 1)] for s < . 
If we use the fact that fj(0) = since f2(0) = 1, and subtract the function P(s) 
from Eq. ([25]), we find 

J[S) J[S) + tt hml ** s'(s> - s)(s> - 4miy/z [Z ° 

which involves the integration of the imaginary part of the function fij 1 (s) over 
the right-hand cut only. This allows us to use hadronic information for s > Am\ 
to constrain the input into our description of the limited 77 — > tttc experimental 
results. 



The factor y 1/3 in Eq. ([22|) is the appropriate Clebsch-Gordan coefficient, as ob- 
tained by decomposing the amplitudes with definite isospin in terms of the ampli- 
tudes with definite charge quantum numbers. As a consequence of the fact that the 
7T°7r Born amplitude is zero, we have: 

5(77 - ttV) = ^|^= 2 - ^-B*=» = , (28) 

from which 

B I=0 = V2B I=2 . (29) 
Applying this to the Born amplitude 

B = 5(77 - vr+TT-) = \/f 5 /= ° + \/f 5 /=2 (30) 

we have 

B-~> = JIb, B'=° = f-B. (31) 



• 7 = 05' and D waves: For these, an identical procedure can be used to calculate 
the modifications from final state interactions reliably up to 0.6 GeV, starting from 

*jf00 = f^Pfx{s)^jf{s) , (32) 
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and using information about the I = S and D waves wir phase-shift 03, 131 to 



compute the f2's. Below 600 MeV, the effect of the unknown 77 — > tttt isoscalar 
phases in the inelastic regime above 1 GeV is small — see Ref. [p3fl . The effect of 
the modification from final state interactions is shown for example in Fig. 3. At 
higher energies, non-pion exchange contributions become increasingly important as 
discussed in || and so these partial waves can no longer be reliably calculated from 
first principles. Instead, we leave the data to determine the amplitudes using the 
representation given by Eq. (7). 

Let us summarize the input and the constraints, in the region under study. Everywhere, 
all the 1 = 2 partial waves and the 1 = waves with J > 4 are given by the modified Born 
amplitudes. Final state interactions only appreciably affect the I = S and D waves. 
The 1 = waves with J = 0, 2 can be reliably predicted by the modified Born amplitudes 
below 600 MeV. However, everywhere they can be represented by Eq. (7) which follows 
merely from coupled channel unitarity. 

With the hadronic amplitudes Ts and 7p known, the four coupling functions af-, a DX 
are what the data and the above constraints determine. These four a's are polynomials 
in energy, which we allow to be at most cubic: this gives a reasonable degree of flexibility, 
but without an overdue number of unphysical structures that could affect the reliability 
of the fit. They are written as a Legendre expansion in terms of the variable x, defined as 

IE — E\ — Ei 

x = F F > ( 33 ) 

so that the energy interval {E\ = 0.28, E 2 = 1-44}, which gives the boundaries of the 
energy range we are fitting, maps onto the interval — 1 < x < 1. Symbolically, we write 
each a as 

3 

a, = 5>i n) P n (x) . (34) 

n=0 

Before we consider the data we are going to analyse to determine their partial wave 
content, let us stress that there is an important region below 600 MeV, where the 77 — ► 
amplitudes are predicted and must also agree with the unitary representation of Eq. (7). 
These predictions provide a reliable starting point for such a general representation. We 
next describe how we build in this constraint. 
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3 Analysis procedure 



3.1 Low energy inputs 

To implement the constraint from Low's theorem |I0[, that as we go down in energy the 



amplitudes are given by their Born terms [11|, we adopt the following strategy. 
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Figure 3: The amplitudes Fjx in the low energy region. The shaded areas rep- 
resent the bands we give as a constraining input, by building a "horn" around the 
central curve generated by the final state modified Born amplitudes. The solid 
lines are the output of our fit. For this plot we have used, as an illustration, the 
outcomes corresponding to what is later called the peak solution; but the difference 
between the two solutions in the low energy range is hardly noticable. 

We take the [.^ja^Iow energy below 600 MeV, as calculated in Sect. 2. Our amplitudes, fitted 
to experiment, must agree with these within some tolerance. To fix this, we construct a 
"horn" around the curves given by by assigning error bars, which are zero at 

threshold and become progressively larger as the energy increases. This reflects the fact 
that neglecting other exchanges then the pion becomes a poorer approximation as the 
energy increases. To do this we introduce a contribution to the total x 2 from fitting 



Xi 



.ov energy 



r W=oi _ \nrl=0] \ 2 

|y J\ J trial form |y J\ \ low energy 1 

AT J 



(35) 



where A J 7 is obviously given by the errors on [^j A °]iow energy we have introduced. These 
ensure an extremely tight constraint close to mi threshold, while allowing larger flexibil- 
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Experiment 


Process 


Int. X-sect. 


| COS 1 m ax 


Ang. distrib. 


| COS 1 max 


Mark II 


77 — > 7T + 7T~ 


87 


0.6 


69 


0.6 


Cr. Ball 


77 — > 7T 7T 


26 


U.o ^OrjooJ 

0.7 (CB92) 


80 


0.8 


CELLO 


77 — > 7T + 7T~ 


30 


0.6 


127 (Harjes) 
249 (Behrend) 


0.55 - 0.8 



Table 1: Number of data in each experiment below 1.4 GeV. Mark II results are from 
Boyer et al. [Rf, Crystal Ball from Marsiske et al. (CB88) M and Bienlein et al. 



(CB92) |3|, and CELLO from Harjes g| and Behrend et al. p 



ity when the two photon energy approaches 500 — 600 MeV. Fig. 3 illustrates how the 
amplitudes T 1 ^ determined by our fit fulfil the low energy constraints given by the Born 
amplitudes modified by final state interactions of Eq. (32). The shaded regions are the 
"horns" inside which the low energy amplitudes are required to fall. 



3.2 Data analysis 

As mentioned in the introduction, the data-sets on two photon scattering into charged pion 
final states at low energies come from Mark II H] and CELLO [T2| , whereas two different 



runs of the Crystal Ball experiment, the last with much higher statistics, provide the 
only available normalized experimental information on 77 — > tt°it for such low energies. 
Table 1 shows the number of data in each experiment, below 1.4 GeV. 

Though the angular distributions contain information about the integrated cross-section, 
because of different bin centres and sizes, these are not always the same. For instance, 
Mark II gives the angular distributions for | cos#| < 0.6 in energy bins of 100 MeV, but 
present the integrated cross-section in 10 MeV steps above 750 MeV. 

From the CELLO experiment, we have angular distributions in Acos# bins of 0.05 from 
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Behrend et al. [I5| and Acos# bins of 0.1 from the thesis of Harjes PH , both in energy 
bins of 50 MeV width. Though these come from the same data sample we believe, we 
have fitted them as separate data-sets but weighted appropriately (see later), since the 
different binning produces quite a difference in the scatter of the data-points, see Fig. 4. 
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Figure 4: Comparison between the angular distribution of CELLO data-points on 
77 — > 7t + tt~ from Behrend |T2j (squares) and Harjes [^4J (diamonds) for some illus- 
trative energy bins (E is in GeV). 



Where statistical and systematic errors are quoted (see tables in 0), these have been 
added in quadrature. Each experiment has an absolute normalization for the cross- 
sections. However, these inevitably provide additional systematic uncertainties. Such 
uncertainties have been included in the results produced by the special low energy trig- 
gering of Mark II. However, above 700 MeV a systematic shift in normalization is apparent 
between the Mark II and CELLO integrated cross sections, though both of them are for 
| cos 9 1 < 0.6, see Fig. 5. It is clear that we must allow for some systematic shift in 
normalization, if we are to describe both data-sets in a sensible way. Mark II quote a 
systematic normalization uncertainty of 7%. With this in mind, we allow for up to a 5% 
relative shift in normalization between Mark II and CELLO experiments. 

For the 7r°7r channel, Crystal Ball had two distinct runs. The first covering the energy 
region from nn threshold, called CB88 [0. The second had 1.5 times as much data, but 



only above 800 MeV. Bienlein et al. |Kj combined this with the CB88 set to produce their 
complete CB92 dataset above 800 MeV. Nevertheless, there was always a clear systematic 
difference in the earlier and later runs through the /2-region. The CB88 set had a higher, 
narrower peak than the combined CB92 set (see Fig. 6). Since the CB88 set (above 800 
MeV) is subsumed in CB92, we cannot really separate them. Consequently, we allow for 
a systematic shift of 3% between CB88 below 800 MeV and CB92 above. 
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Figure 5: Comparison between Mark II [§] (diamonds) and CELLO [|T^] (squares) 
integrated cross-sections for 77 — > 7r + 7r~ in their common energy range. 
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Figure 6: Comparison between CB88 [0 (diamonds) and CB92 |l3j (squares) 
integrated cross-sections for 77 — > 7r°7r°. 
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While both charged and neutral pion experiments would be consistent (and the fits even 
better) with larger normalization shifts, we have erred on the side of caution in allowing 
this freedom. To repeat : we assume that the results from CELLO, CB88 below 800 MeV 
and Mark II below 400 MeV are absolute, and allow small shifts of other data-sets with 
respect to these. 

3.3 Weighting, relative normalizations and fitting 

From Table 1, we see that the number of data-points is far greater for the reaction 
77 — > 7r + 7r~ than for 77 — > 7r°7r°. Since an accurate separation of the 1 = component 
requires both are accurately described, we give different weight factors to each dataset. 
We choose these so that the Mark II and the CELLO data have roughly the same number 
of weighted data, while the weight assigned to the Crystal Ball data approximately equals 
the weighted sum of the Mark II and CELLO data. Nevertheless, good agreement is not 
easy to achieve. 

The analysis program (GAMP ||) works by integrating the amplitudes over the appro- 
priate bin in energy and angle for each data-point. It does not just use the central values. 
This is to allow for any strong local variation of the amplitudes, particularly near KK 
threshold. In the next Section, where we display the solutions we find, this should be 
borne in mind. Where the energy bins are sizable (as with Crystal Ball), histograms are 
plotted (see Figs. 8 and 11). Where the energy bins are fine (as with Mark II data in 
10 MeV steps), the fits are shown more appropriately as continuous lines joining the bin 
centres, as in Figs. 7, 9, 10, 12. However throughout, the fits are histograms. 
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4 Results 



Our fits deliver two classes of distinct solutions. As is seen from Figs. 7-18, the two 
classes have very similar quality as far as fitting the data is concerned, yet they have 
quite distinct characteristics. The first has a peak in the cross-section located in the 1 
GeV energy region and from now on will be referred to as the peak solution. The second 
has a dip in the same energy region and will be called the dip solution. The plots with 
the fits to the Mark II it + it~ integrated cross-section look more structured than those 
of CELLO in the 1 GeV region (cf. Figs 7 and 9, or 10 and 12). This is because our 
fitting routine integrates over bins, which for Mark II are only 10 MeV wide in this region 
compared with 25 MeV from CELLO. Detailed dynamical features are picked up more 
strongly in the finer binned data. Recall the fits are not continuous curves but histograms, 
as shown in Figs. 8, 11 for the -K n° data of Crystal Ball. 

Looking at the plots of just the Mark II results, Figs 7, 10, on the integrated cross-section 
for 77 — > 7r + 7r~ in the 1 GeV region, one might be tempted to conclude that the peak 
solution is disfavoured. However, one cannot conclude that from the CELLO results, 
Figs. 9, 12, on the same channel and Crystal Ball, Figs. 8, 11, on the tc°tt° final state, nor 
indeed from looking at the fits over the whole energy region for Mark II. Thus, individual 
features are a poor guide, even though one's eye naturally picks those out. Indeed, the 
overall quality of the fits in each sector for the two distinct solutions are quite comparable 
as can be seen from Tables 2 and 3. There we report the contributions to the total x 2 
from each experimental set, for the integrated cross-section and the angular distribution 
separately for the best of these solutions. This total is expressed as the x 2 P er degree of 
freedom. The peak solution has a slightly lower overall x 2 ■ However, the CELLO data 
turn out to be easiest to fit (and so we have reduced their weight by 1/2 to achieve greater 
sensitivity to the rest). In contrast, the two solutions have appreciably greater x 2 for the 
Mark II and Crystal Ball results by ~ 0.5 per degree of freedom. 

To explore the neighbourhood of these best fits, it is convenient to characterise these 
classes of solutions by the relative amount of / = S, D and D 2 contribution to the 
cross-section a = a s + o"d„ + o"d 2 at the / 2 (1270) peak. Each solution corresponds to a 
point in the equilateral triangle of height 1 with sides os/cr = 0, o^fo = and o\d 2 /o" — 
(see Figs. 15, 16). For the two classes of solutions, peak and dip, we display the overall x 2 
per degree of freedom found for each fit in these equilateral triangles. 
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PEAK SOLUTION X L = 1-40 


DAUCi llllCLlu 






Y 2 

A average 




A n cr Hit;tT*iVi 


Mark II 


77 — > 7r + 7r~ 


156 


1.54 


X 2 = 1-82 


X 2 = 1.19 


Cr. Ball 


77 — > 7r°7r° 


106 


1.44 


X 2 = 1-42 


X 2 = 1-44 


CELLO 


77 — > 7T + 7T~ 


406 


1.33 


X 2 = 0.65 


X 2 =1-13 
from Harjes 
X 2 =1-52 
from Behrend 



Table 2: Summary of contributions from each experiment to the total v 2 for the peak 
solution. Here Xtot i s calculated by dividing the sum of the x 2 's for all data-sets by 
the total number of data-points we are fitting, namely 668. Xaveragc i s computed by 
dividing the sum of the x 2 's for each data-set by the total number of data-points in 
that experiment. 



DIP SOLUTION xL = !-48 


Experiment 


Process 


data-points 


Xaverage 


Int. X-sect. 


Ang. distrib. 


Mark II 


77 — > 7T + 7T~ 


156 


1.75 


X 2 = 1-99 


X 2 = 1.46 


Cr. Ball 


77 — > 7r°7r° 


106 


1.62 


X 2 = 1-97 


X 2 = 1.51 


CELLO 


77 — > 7T + 7T~ 


406 


1.33 


X 2 = 0.88 


X 2 =1.09 
from Harjes 
X 2 =1.51 
from Behrend 



Table 3: Summary of contributions from each experiment to the total v 2 for 
the dip solution, as described for Table 2. 
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Figure 7: Cross-section as a function of the tttt invariant mass integrated over 
| cos 9 | < 0.6, and angular distributions as a function of cos 9, for the 
77 — > 7T + 7r _ process, from the MARK II experiment || (peak solution). 
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Figure 8: Cross-section as a function of the tttt invariant mass integrated over 
| cos 9 | < 0.8 for E < 0.8 GeV and | cos 9 \ < 0.7 for higher energies, and 
angular distributions as a function of cos 9, for the 77 — > 7r°7r° process, from 
the CRYSTAL BALL experiment 0, [L3| (peak solution). 
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Figure 9: Cross-section as a function of the tttt invariant mass integrated over 
| cos 9 | < 0.6, and angular distributions as a function of cos 9, for the 77 — > 7r + 7r~ 
process, from the CELLO experiment (peak solution). The angular distributions are 
from the binning of Harjes p4] - see Fig. 13 for the Behrend et al. binning. 
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Figure 10: Cross-section as a function of the irn invariant mass integrated over 
| cos 9 | < 0.6, and angular distributions as a function of cos 9, for the 
77 — > 7T + 7r _ process, from the MARK II experiment || (dip solution). 
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Figure 11: Cross-section as a function of the tttt invariant mass integrated over 
| cos 6 | < 0.8 for E < 0.8 GeV and | cos 9 \ < 0.7 for higher energies, and 
angular distributions as a function of cos 9, for the 77 — > 7r°7r° process, from 
the CRYSTAL BALL experiment 0, [I3j (dip solution). 
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Figure 12: Cross-section as a function of the nir invariant mass integrated over 
| cos 9 | < 0.6, and angular distributions as a function of cos 9, for the 77 — > 7r + 7r~ 
process, from the CELLO experiment (dip solution). The angular distributions are 
from the binning of Harjes p4] - see Fig. 14 for the Behrend et al. binning. 
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Figure 13: The angular distributions as a function of cos 6, for the 77 — > 7r + 7r 



process, in different 7T7t mass bins from the Behrend et al. analysis [12] of 
CELLO experiment for the peak solution, to be compared with Figs. 9, 12, 
where the bins in cos# are larger. 
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Figure 14: As for Fig. 13, but showing the dip solution. 
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These diagrams, Figs. 15, 16, clearly show that the fit singles out a well defined region 
in the parameter space for each class. If we tried to drive them outside this area, their 
X 2 values would increase very rapidly. Each rectangular flag in Figs. 15, 16 indicates 
one solution is found at that particular point, and the number which labels it is the 
corresponding \ 2 . The round labels indicate the position of solutions A (the typical 
"good" solution) and B (technically the "best" solution) found by Morgan and Pennington 
in their previous 77 — > tttt data analysis |§. Notice how, while their favoured solution A, 
falls inside the region determined by both our classes of solutions, their best solution B, 
having very large S and D wave contributions at the /2(1270) peak, lies considerably far 
away from it. 

If we now compare the diagrams in Figs. 15 and 16, we see that our peak class of solutions 
singles out a region in the parameter space where the as/a components are remarkably 
small, approximately bounded by 



Alternatively, the dip class of solutions determine a region the boundaries of which are 
given by higher values of as /a and lower values of <7d /(x: 

0.11 < a s /a < 0.35, 



Furthermore, in Tables 4 and 5 we report the x 2 ' s of two representative solutions from 
the peak and dip classes, to be compared with those corresponding to our two favoured 
solutions, respectively, shown in Tables 2 and 3. They have been picked out from the 
solutions lying in the regions of the parameter space shown in Figs. 15 and 16 in order to 
show the difference in the x 2 's corresponding to each data set for an ordinary solution in 
either the peak or the dip class compared to our best solution (indicated by the shaded 
flag). It is immediately evident how the x 2 ' s corresponding to the CELLO experimental 
data, not only the overall averaged one but also the individual x 2 for the integrated cross 
section and the angular distribution, hardly change at all. On the contrary, big variations 
occur for the other data-sets, the x 2 °f which in some cases increase by a factor of 2. 
Indeed, it is all the data-sets together that determine the features of the solutions, but 
we once again want to stress the point that while a remarkably good agreement with 
the CELLO data is easily achieved most of the time, the Crystal Ball and Mark II data 
are always hard to satisfy simultaneously. In fact, as seen from Figs. 9 and 12, the 



0.00 < a s /a < 0.25 , 
0.06 < a D Ja < 0.35 . 



(36) 



0.03 < a D Ja < 0.16 . 



(37) 
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Figure 16: Mapping of an assortment of solutions with different proportions 
of as, &d and <td 2 at 1270 MeV labelled by the x 2 per dof for the dip class. 
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CELLO data have finer energy and angular bins. Consequently, it is these data that most 
powerfully constrain our solutions. Thus the solutions in Tables 2-5 all have very similar 
X 2 for this sector, and it is in the contributions to x 2 from the Mark II and Crystal Ball 
data-sets that they differ. A final comment ought to be made on the contribution to the 
total x 2 from the low energy region constraints. As we have seen in Fig. 3, the agreement 
of the fit solutions to the low energy amplitudes Fjx , calculated by a dispersion relation 
over the right-hand cut, is pretty good. The ratio of the low energy x 2 , see Eq. fl35|), to 
the total overall x 2 is the following 

XL energy /Xtot = °-22 for the peak solution , (38) 
XL energy /xL = 0.17 for the dip solution . (39) 



Illustrative solution no. 1 Xtot = 1-82 


Experiment 


V 2 

A average 


Int. X-sect. 


Ang. distr. 


Mark II 


2.82 


3.21 


2.32 


Cr. Ball 


2.30 


3.09 


2.03 


CELLO 


1.30 


0.64 


1.35 



Table 4: Summary of contributions from each experiment to the total x 2 f° r 
an illustrative solution in the peak class (see Fig. 15) to be compared with the 
best solution in that class. See text for comments. 



Illustrative solution no. 2 Xtot = 2-00 


Experiment 


V 2 

A average 


Int. X-sect. 


Ang. distr. 


Mark II 


2.87 


2.87 


2.87 


Cr. Ball 


2.90 


3.71 


2.64 


CELLO 


1.42 


0.63 


1.49 



Table 5: Summary of contributions from each experiment to the total x 2 f° r 
an illustrative solution in the dip class (see Fig. 16) to be compared with the 
best solution in that class. See text for comments. 
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Figure 17: / = 77 — > 7r7r cross-section as a function of energy (peak solution). 
Notice the two pronounced peaks in the S partial wave: the first just above thresh- 
old, where it dominates the 7r + 7r~ cross-section, and the second corresponding to 
the / (980) resonance. 
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Figure 18: / = 77 — > irir cross-section as a function of energy (dip solution). 
Notice the pronounced peak in the S partial wave at threshold, where it dominates 
the 7r + 7r~ cross-section, and the dip corresponding to the / (980) resonance. 
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5 Solutions 



Figs. 17 and 18 show the dominant partial wave components of the I = 77 — > nn cross- 
section for our two favoured solutions. These fall into two classes characterised by having 
either a peak or a dip in the S'-wave cross-section in the 1 GeV region, corresponding 
to two distinct coupling patterns for the /q(980). In the peak solution Ort(s) is larger 
than a n (s) in the 1 GeV region so that, in the decomposition given by Eq. (7), the 
contribution of the hadronic reaction tttt — > KK dominates over titt — > titt (see Figs. 17 
and 2). As a consequence the /o(980) has a larger coupling to 77 — > KK in the peak 
solution. However, there is nevertheless a crucial contribution from a n (s), and hence of 
the 7T7T — ► 7T7T amplitude, which results in the /o(980) peak in this 77 amplitude (Fig. 17) 
being sharper than that seen in titt — > KK (Fig. 2). In the alternative class of solutions, 
o£(s) dominates over a|(s) in the 1 GeV region and as is clear from Fig. 18 the /o(980) 
then manifest itself as a dip as in the titt — > ttti cross-section (see Fig. 2). 

In the energy region above 1 GeV, the cross section is dominated by the D 2 -wave embodied 
in the /2(1270) resonance. An oversimplified description of individual channels is to 
ascribe the peak in this region wholly to /2(1270) formation in the A = 2 state |J. Such a 
simplification is often used when trying to extract a 77 width for the tensor mesons from 
a single charge final state with limited angular coverage, when a true amplitude analysis 
is not possible. Here, as in the earlier analysis by Morgan and Pennington |J, sizable 
contributions of both S and Dq waves in this region are strongly favoured. This is in 
good agreement with solution A of Ref. ||: however we do not find such large S and 
-D contribution as their "technically best" solution, B. We report here Os/cr and Udq/ 17 
ratios for the 1 = cross-sections at 1270 MeV for the best solutions in each class: 



o s ja = 0.13 
<t d J<t = 0.22 

a s /a = 0.25 
&D /a = 0.13 



> for the peak solution , (40) 



> for the dip solution . (41) 



Notice in Figs. 17, 18 how the influence of the Born term at low energies (modified of 
course by crucial final state interactions) means that even outside the /o(980) region the 
S'-wave cross-section is not simply describable by a Breit-Wigner resonance at 600-1200 
MeV, rather its contribution is spread over a wide region. 
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5.1 77 couplings 



We now calculate the 77 couplings of the resonant states in the threshold to 1.4 GeV 
region that our amplitude analysis determines. We do this by two different methods ||. 
The first is based on the analytic continuation of the amplitudes we have found in our 
fit into the complex s-plane to the pole position. This is the only formally correct way 
of deducing the couplings of any resonance and its outcomes are free from background 
contamination. The second is a more naive approach based on the Breit-Wigner-like peak 
height. For the / 2 (1270), these two methods give nearly identical results, as expected for 
an uncomplicated and isolated resonance with a relatively nearby second sheet pole. For 
the /q(980) only the pole method is applicable because of the overlapping of this state 
with KK threshold and the broader /o(400 — 1200). For this latter broad state, only 
the "peak height" provides a sensible measure of its 77 width, since its pole is too far 
from the real axis to be reliably located under the approximations needed to perform the 
analytic continuation, as will be clear from what follows. 

To work out the pole residue based definition of the radiative widths, we suppose the 
strong interaction amplitudes Tj(s) and the corresponding 77 — > amplitudes J-'(s) to 
be dominated by their pole contribution near the resonance pole, i.e. for s ~ sr; then we 
can write them in the form 

q 2 

Tj(nn -> 7T7r)(s ~ s R ) = — ^— , (42) 

sr — s 

-> i™){s ~ s R ) = 9l9lz . (43) 
sr - s 

It is easy to see that the couplings g n and g 7 can be extracted as the residues of these 
amplitudes at the resonance pole sr. Thanks to the parametrization of the amplitudes T 
in terms of the K— matrix elements, as given in Eq. (10), we know their numerator, N„(s) 
and Nk{s) respectively, and denominator, D(s), which is the same for the two of them. 
So we can use the expressions 

Tjiltlt — > 7T7r) 

Tj(nir -> KK) 

and Eq. (7), to write J 7 ja(77 — » vr7r) as: 

aj\s) N n (s)+a J K \s)N K (s) 
*>a(77 - tttt) = — 



D(s) 

N K (s) 
D(s) 



(44) 
(45) 
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Now, in the region nearby s R we can make a Taylor expansion of the function D(s), 
truncated at the first order 



D(s ~ s R ) ~ D(s R ) + D'(s R )(s - Sji ) = D'(s R )(s - s R ) , (47) 

where by definition D(s R ) = at the resonance pole. Finally, by substituting this into 
Eqs. (|44]J45|) and comparing with Eqs. (f42|J[3"l) , we find 



2 K(S R ) 



r? 7 ( P ole) = , (50) 



* = ^ ■ <«> 
and 

from which the coupling g 1 can readily be calculated. The corresponding width is then 
evaluated using the formula 

® 2 Pr |ff 7 | 2 
A(2J+l)m R 

where a = 1/137 is the fine structure constant, and /3 R = (1 — 4m^/m|,) 1 / 2 ~ 1. 

Because of its very large width, the /o(400 — 1200) coupling to two photons cannot be 
calculated with this technique, since one cannot reliably continue so far into the complex 
plane. As an alternative, we give a rough estimate of its width by using an expression 
based on the standard resonance peak formula 

Pit^u _ craves, peak) m 2 R r tat 

where BR is the hadronic branching ratio for the final state considered. 

Table 6 shows the results we obtain for the 77 widths of the / (980), /o(400 — 1200), 
/2(1270) from either the pole or the peak-height definition, as appropriate, when choosing 
solution 1 and solution 2, respectively. 

The continuation to the second sheet pole for the /o(980) is rather sensitive to the 
parametrization of the i^-matrix (and hence of the T-matrix, Eq. 10) used. In the 
Re VAMP analysis Jl9| described in Sect. 2, the i^-matrix elements are given by a single 



pole plus different order polynomials. These each fit the hadronic data equally well. For 
the dip solution, the 77 width for the / (980) is 0.31 keV for ReVAMPl, and 0.34 keV 
for ReVAMP2, so rather little change. However, for the peak solution for which there is a 
much stronger interplay between the 7s(7T7r — > 7T7t) and the Ts(7rn — > KK) contributions 
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r(/ 2 (1270)) 


r(/ (980)) 


r(/ (400 - 1200)) 


Peak solution 


3.04 


0.13-0.36 


3.0 


Dip solution 


2.64 


0.32 


4.7 



Table 6: Two photon partial widths in keV for the meson states in 
our two classes of solutions. 



to Eq. (7), we find a 77 width of 0.13 keV for Re VAMP 1 and 0.36 keV for ReVAMP2. 
Hence the values in Table 6. 

From the variation between different solutions within the classes indicated in Figs. 15 
and 16, we estimate that the uncertainty on the 77 width of the ^(1270) is ±0.08 keV 
within each class of solutions. To this must be added a 5% uncertainty in the absolute 
normalization. Consequently, it is the difference between solutions (rather than within a 
given class) that constitutes the major uncertainty, and so we quote 

r(/ 2 (1270) -> 77) = (2.84 ± 0.35) keV . (52) 

Whilst this is in good agreement with the PDG'98 estimated value of (2.8 ± 0.4) keV, 
it is somewhat at variance with the PDG'98 fitted value of (2 

.44ta29j keV, based on 

several different analyses M of either 7c + 7T~ or 7r°7r° data separately using quite different 
assumptions. It is important to stress that our value is nearest to a model independent 
amplitude analysis result one can presently achieve. 

The uncertanties on the 77 widths of the /o(980) and /o(400 — 1200) are more problematic. 
For the /o(980) we quote 

r(/ (980) - 77) = (0.28 keV. (53) 

However, a decision on whether the dip or peak solution was correct would reduce the 
uncertainty dramatically. 

For the /o(400 — 1200) the estimate is much cruder and a 50% uncertainty is likely, giving 

r(/o(400 - 1200) -> 77) = (3.8 ± 1.5) keV. (54) 
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Once again discriminating between the classes of solutions would change the central value 
and reduce the error within these ranges. 

What about the composition of these states ? Quite independently, we, of course, know 
the /2(1270) is the rin member of an ideally mixed multiplet. The 77 widths of the neutral 
members are then expected to be in the ratio of the squares of the average squared charges 
of their constituents, so 

r(/ 2 -> 77) : T(a 2 - 77) : T{f 2 - 77) = 25 : 9 : 2 . (55) 

Experiments give ratios very close to this. The prediction for the 77 width not only 
depends on the charges of the constituents to which the photons couple, but to the 
probability that these constituents annihilate — their overlap. For members of a qq 
multiplet we expect these probabilities should be roughly equal. Experiment for the 
lowest tensor nonet confirms this. (Any differences can readily be explained by a small 
departure from ideal mixing — see the Appendix of Ref. |§.) 

For the scalars, we need a certain modelling. The simplest would be to assume the lightest 
scalars are the shadow of the tensor nonet. Then non-relativistically, Chanowitz [p5|| 
deduced that the corresponding tensor and scalar two photon widths are related by 

15\ / m \ 3 



r(o++ - 77) = [f) x j x r(2++ -> 77) , (56) 

where the factor of (mo/7712) 3 takes account of the mass splitting. This gives the pre- 
dictions in Table 7, under the assumption that the / states are either purely ss or 
non-strange nn. Relativistic corrections to Eq. (56) have been computed [^6j and found 
to reduce the ratio r(0 ++ — > 77) /T(2 ++ — > 77) by as much as a factor of 2 for light quark 
systems. An oft-cited alternative structure for the /o(980) IS cLS du KK molecule ||15|| . 



Then though the fourth power of the charge of the constituents is far greater for a KK 
molecule than for a simple ss bound state, the molecule is a much more diffuse system, 
so that the probability of the kaons annihilating to form photons is strongly suppressed. 



Thus in the non-relativistic potential model of Weinstein and Isgur ||15|| , computation of 
the molecular radiative width gives F(f (KK — > 77)) = 0.6 keV P7| . 



Even this may be far too simplistic for these ++ states. From the work of Tornqvist Jjj 



(replicated by one of us [28]), we know the dressing of bare qq states, by the hadrons into 
which the physically observed mesons decay, is particularly large for scalars. In lattice- 
speak, their unquenching is a big effect. In a scheme where it's assumed pseudoscalar 
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nn 


ss 


KK 


r ++ 


4.5 


0.4 


0.6 



Table 7: Two photon partial widths in keV predicted for a 
conventional qq nonet of scalar states and for a KK molecule. 



meson pairs provide the dominant dressing, the /o(980) has not only an ss (and smaller 
nn) component, but a large KK admixture, too. It is not a case of the physical hadron 
being either a qq bound state or a KK molecule, but it is in fact both! What this non- 
perturbative treatment predicts for the radiative width of the /q(980) is a calculation 
under way. 

We know there are other scalars beyond the region of our analysis: the /o(1500) and 
/j(1710) in particular [Q. Some claim there is also an /o(1370). In the past, we have 
argued that this may just be the higher end of the broad /o(400 — 1200), cut-off below 
by KK and An phase-space. Interestingly, the very recent potential model analysis of 
Kaminski et al. |29[ suggests the / (1370) may actually be the same as the / (1500). In any 
event, the radiative width of all these states is a key pointer to their composition, gluonic 
or otherwise [17, 30|. At present, only crude estimates are possible, for instance [3~I] that 
r(/o(1500) — > 77) < 0.17 keV. To achieve something more is the challenge for the future. 
A major task is to extend the present Amplitude Analysis beyond 1.4 GeV. This requires 
a study of two photon production of not just two pion final states, but 4ir and KK too. 
Only by a detailed analysis of these final states simultaneously, can we hope to extract a 
true scalar signal from under the dominant spin 2 effects in the region from 1.3-1.8 GeV 
and so deduce the radiative widths of the / (1500) and /j(1710). 

We have seen that presently published data allow two classes of solutions, distinguished by 
the way the /o(980) couples to 77. A primary aim must be to distinguish between these. 
We believe that data with sufficient precision may already have been taken at CLEO that 
could do this p2| . However, these results have not yet been corrected for acceptance and 
efficiency, and sadly may not be. To go to higher masses within the resonance region may 
well be possible when corrected results from LEP2 and future experiments at £>-factories 
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become available. The challenge for theory is equally demanding: it is to deduce what 
the radiative widths for the /o(980) and for the /o(400 — 1200) we have determined here 
from experiment, and summarised in Table 6, tell us about the underlying nature of these 
dressed hadrons. Only then may we hope to solve the enigma of the scalars: states that 
are intimately related to the breakdown of chiral symmetry and hence reflect the very 
nature of the QCD vacuum. 
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